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Abstract 

The phonon dispersion of chiral single wall carbon nanotubes has been obtained from 6x6 
dyanimic matrix. The present manuscript is the extension of G. D. Mahan and Gun Sang Jeon's 
work on armchair and zigzag nanotubes. [see P. R. B., 70, 075405 (2004)]. We use spring and 
mass model with the proper phonon potential suggested by ^3). We can calculate the dispersion 
of single wall carbon nanotubes near T point with arbitrary chirality. The results are compatible 
with Mahan et. al's results for armchair and zigzag SWNTs. 
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Carbon nanotubes have attracted much interest on their electronic and thermal transport 
properties [l]. Single wall carbon nanotubes (SWNTs) can be thought of as a sheet of 
raphite rolled into a cylinder Q. In many experiments, phonons play an important role 
a, |3j , especially, many of properties of SWNTs are determintd by the low-energy phonon 
excitations. Therefore there is a pressing need for understanding the phonon dispersion in 
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)on nanotubes, There are several works on this topic 
"1 

20|| . G. D. Mahan and Gun Sang Jeon proposed a tractable method in terms of spring 
and mass model to calculate the phonon dispersion of armchair and zigzag SWNTs [17 1. 

In the present paper, we investigate theoretically the phonon dispersion of SWNTs with 
arbitrary chiral index (n,m). We follow the algorithm proposed in [l7| and we extend it for 
our calculations. 

In SWNT, there are two carbon atoms in a unit cell. We label A atoms at the lattice 
site Rj and their three nearest neighbors B atoms at site + 6. Unit vector 8 reflects the 
direction of S. 

FigG] shows the commonly used coordinate system in nanotubes with base vectors 

a x = V3a(V3/2, 1/2), 
a 2 = ^(^3/2, -1/2),, 

The chirality of a certain carbon nanotube is characterized by a pair of integer (n,m). 
When a sheet of graphite is rolled into a cylinder, the points (0,0) and (n,m) are identical, 
and the diameter of tube R is the distance between above two points on the 2D plane divided 
by 2tt. 

Consider a type A Atom located at lattice site (t,l) on 2D graphite sheet. For an 
(n,m) SWNT with radius R, the atom position in 3D Descartes coordinate system is, 
(c = 3o/2, 9 Z = V3a/R) 

R A = [Rcos(9 u ),Rsm(9 tl )Jc], 
R% = [Rcos(6 u + Ui8 z ),Rsm(8u +Ui9 z ),lc + avi], 
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i = 3. 
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Where I = & ■ . l 2 n ~ tm .. 



2 V n 2 +nm+m 2 ' 

denotes the position along the z axis. 



It is easy to verify that the carbon atoms' coordinates of armchair tube m = n and zigzag 



tube m = shown in ref. 
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171 ] are just special examples of our general expression. 



According to the atom coordinates given above, the three unit vectors Si are 



#0 = [-V3 Ul sm(8 tl + %9 Z ), V3 Ui cos(0 t i + ^9 
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1,2,3. 
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We know that for carbon nanotubes, the basis of lattice vectors in real space, in cylindical 
coordinate system, 
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a 2 = (0,T). 
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The corresponding reciprocal lattice vectors are 



Here, 



It is easy to verify that, 



/3i = (iV,0), 
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<ifj = gcd(n,m). 



OLi ■ (3 a = 2ir5ij. 



In SWNTs, the electron states have two quantum numbers (q, a), where continuous index 
q is the wave vector along the tube axis. Discrete quantum number a reflects the angular 



dependence around the tube. In first Brillouin zone, then, 

a = 0,l,...,iV-l, -l< q <l. 

In cylinder coordinate system [p, 9, z], we consider the displacement of the atoms A and 
B in three directors Q sp , Q s g, Q sz , (s = A, B) are, 

Qa,u = Q A ,ti,p[ cos ^ti),sm(9 t i),0} 

+ Qa,uA- Bin(0tf), cos(0 ti ), 0] + Q A>U , Z [0, 0, 1], (6) 

Q%i = QB,ti, P [cos{9 t i + Ui6 z ), sin(0« + Ui9 z ), 0] 

+ Qb,u,o[- sia (9 t i + Ui6 z , cos(9 u + Ui9 z ), 0] 

+ Q B , tl , z [0,0,1], 1 = 1,2,3. (7) 

The wave vectors occur in the collective coordinates. Wave vector q is used for the tube 
axis, the other quantum number a denotes the position around the circumference of the 
tube. Similar to the discussion for the single electron dispersion, we have a = 0, 1..., iV — 1 
with N = 2( " 2+ " m+m2) - Thus, 

Qa,u,t = —7= QA, qa ,r exp[i(qcl) + a9 u ], 

V-'*a qa 

Q { B,ti,r = -7== Q b ,qa,r exp[i(qd + iqaVi) + a(9 u + Ui9 z )), 

V-'*a qa 

where r = (p, 9, z) and N a is the total number of carbon atoms in SWNTs. 
The phonon potential energy Vi is [l7| 



vi = ^£EI£ (i) i 2 , (S) 



A 

= e i te d+pBtl) E®,i = 1,2,3 (9) 
S(<) = e ^^».(V3 Ui sin(^)Qfl >p + V3u i cos(^)Q B e + «igB*) 

+ v / 3uiSin(^^)Q j4/9 - v / 3n i cos(^y i )(5A0 - ^Qa 2 , i= 1,2,3. (10) 

The displacement Q is a complex quantity, and the real and imaginary displacements are 
independent variables. According to ref. [l7| . we only need the force on Q as 5V/5Q* to 
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obtain corresponding determinant equation, 



= det 
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where Si Z = sm(^^) , c i2 = cos( 2 ^) and /3, = qav j + aui9 z . For convenience, a dimension- 
less quantity fti 2 = rrvJ 1 jK\ is defined, i^i is the spring constant. 

If the SWNT is zigzag or armchair type, our results are reduced to the results shown in 
ref. The above determinant equation has analytical solutions when q = 0. Since a has 
iV possible values, the phonon modes of chiral SWNTs' number is much more than achiral 

;hree zero-frequency modes 



171] 's discussion, we have to 



SWNTs'. The calculated results are unphysical, for there are 
and have no acoustical modes or flexure modes. Following ref. 
consider the second neighbor interaction. 

We thus use springs between second neighbors. For A atoms with coordinate Ra, its six 
second neighbors are with coordinates 



0(1,2,3) 

It A 



[Rcos(9 u + Vi9 z ), Rsm(9u + Vi9 z ), Ic + 3aw.j], i — 1, 2, 3 
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and the unit vectors to them are 
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We get 
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+ Vi sm(^^)Q Ap + Vi cos(^^)Q Ag + V3~UiQ Az , i = 1, 2, 3. 



(15) 



We reduce above equations to 



2e«/ 2 {^ sin(^)Q A , p cos A/2 + ifo cos(^)Q A , 



+ v^Q^] sin(A/2)} 
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where /3j = 3qaUi + a6 z Vi and the z* denotes the complex conjugate of z. Obviously, 

~(*+3) 

H | = H |. These interactions give a dynamical matrix for the oscillation of the A 



sublattice 
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where, = sin( : ^),c^, = cos( :H y L ), = sin(^) and Cj = cos(y). The dimensionless 
quantity f2 2 2 = Mw 2 / K%. Ki is the spring constant. 

An identical set of equations applies to the second neighbor interactions among the B 
atoms. From the above determinant equation, we find there are two acoustical modes 
at a = 0, i.e., longitudinal acoustical phonon and the twist mode. From corresponding 

n 

velocities of the sound Vl and vtw-, we ma Y fit the spring constants 1]. In addition, when 



a = 1, there occurs one flexure mode. It implies that we should consider radial bond bending 
effect to get ZO phonon modes, which is similar to the out-of-plane mode of graphene. 

For a sheet of graphene, the vibrations perpendicular to the plane are well described by 
bond-bending forces. For SWNTs, in order to investigate the radial bond bending, ref. 
gives the relative phonon potential energy 
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where are the first neighbors of Qj. If Qj denotes an A atom, the normal vectors are 



fii = [cos{8 tl + -£0z), sin(0« + y6» 2 ), 0], 
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In polar coordinates, using the Fourier transformations, these two interactions are 
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The interaction leads to a dynamical matrix 
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where c 



ELcos^ 
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£f =1 sin(^) and s 



E?=i sm(^j^)e ll3i . The dimensionless quantity £l 2 = Muj 2 /2K%. K3 is the spring constant. 

When a = 0, the solution of above equation does have ZO modes. At q — 0, u>zo is also 
important to fit spring constants. 

We add together all three interactions. With the weights r\ = 1, r-z = K2/K1 and 
r 3 = K3/K1, we get six equations: 

V 2 Qap = [(j23u 2 s 2 z )Q Ap + ...] + r 2 [(8j2v^ 2 C 2 )Q Ap + ...] 

8=1 1=1 

+ r 3 [(c 2 + \c\ 2 )Q Ap + •••], etc. 



From these equations, we can calculate the correct phonon dis per sion of SWNTs. To 
varify our method, we calculate a zigzag SWNTs (10,0). As ref. |l2j, we set ft = y/3 to 
the optical phonon at 1600 cm -1 . For this tube, the two acoustic modes have velocities 
vl = 16.4 km/s and vtw — 8.7 km/s. The optical phonons are at frequencies 1600 cm -1 , 
1580 cm" 1 , 877 cm" 1 and 291 cm" 1 . 

In summary, we consider the phonons in chiral SWNTs. Lower symmetry of chiral 
SWNTs make the computation more difficult. We have developed a computational method 
to deal with this difficulty and obtain the reasonable results. It can be used to calculate 
the phonon dispersion of arbitrary chiral SWNTs. From the phonon dispersion, we may do 
further research on the phonons involved phenomena. 
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FIG. 1: The 2D sketch map of single wall carbon nanotubes. This image was made with VMD 
and is owned by the Theoretical and Computational Biophysics Group, an NIH Resource for 
Macromolecular Modeling and Bioinformatics, at the Beckman Institute, University of Illinois at 
Urbana- Champaign . 
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FIG. 2: Phonons in a zigzag (10,0) tube, a = 0,1, 2, N - 1. 
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FIG. 3: Large scale presentation of the region at low frequency and small wave vector of (10,0) 
tube. 
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